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This work has attempted to reconcile puzzling neutrino oscillation results from the LSND, KAR- 
MEN and MiniBooNE experiments. We show that the LSND evidence for z7 M — >• z7 e oscillations, 
its long-standing disagreement with the results from KARMEN, and the anomalous event excess 
observed by MiniBooNE in and data could all be explained by the existence of a heavy 
sterile neutrino (vh)- All these results are found to be consistent with each other, assuming that 
the vh is created in f M neutral-current interactions and decays radiatively into a photon and a 
light neutrino. Assuming the Vh is produced through mixing with v^, the combined analysis of the 
LSND and MiniBooNe excess events suggests that the vh mass is in the range from 40 to 80 MeV, 
the mixing strength is |t/ M h| 2 — 10~ 3 — 10~ 2 , and the lifetime is r„ h < 10 -9 s. Surprisingly, this 
LSND-MiniBooNE parameter window is found to be unconstrained by the results from the most 
sensitive experiments. We set new limits on | C/ M h. | 2 for the favorable mass region from the precision 
measurements of the Michel spectrum by the TWIST experiment. The results obtained provide a 
strong motivation for a sensitive search for the in a near future K decay or neutrino experiments, 
which fit well in the existing and planned experimental programs at CERN or FNAL. The question 
of whether the heavy neutrino is a Dirac or Majorana particle is briefly discussed. 
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PACS numbers: 14.80.-j, 12.60.-i, 13.20.Cz, 13.35.Hb 



I. INTRODUCTION 

Over the past 10 years the LSND collaboration [l| has 
observed an event excess with a significance of 3.8 a at 
LANSCE 0, [|. This excess, originally interpreted as 
a signal from v ^ — > v e oscillations was not confirmed 
by further measurements by a similar experiment KAR- 
MEN which was running at the ISIS neutron spal- 
lation facility of the RAL [Hf. The MiniBooNE experi- 
ment at FNAL 11, designed to examine the LSND effect, 
did not find evidence for —¥ v e oscillations. However, 
an anomalous excess of low energy electron-like events in 
charge-current quasielastic (CCQE) neutrino events over 
the expected standard neutrino interactions has been ob- 
served Q ■ This MiniBooNE anomaly has been confirmed 
by the finding of more excess events [8J. Recently, the 
MiniBooNE experiment has reported new results from a 
search for — » 77 e oscillations Q. An excess of events 
was observed which have a small probability of being 
identified as background-only events. The data are found 
to be consistent with 77^ — » 77 e oscillations in the 0.1 eV 2 
range and with the evidence for antineutrino oscillations 
from the LSND experiment. 

The new observations bring more confusion than clar- 
ity to the experimental situation. The inconsistency be- 
tween the results of the experiments, in particular be- 
tween the LSND and KARMEN experiments, is also 
confusing in light of the apparent simplicity of the pri- 
mary reaction, p(v e ,e)n, used by these experiments for 
the oscillation signal search, and also in view of the fact 
that other results, e.g. the inclusive cross section for 
12 C(v e: e) 12 N* with an electron in the final state, mea- 
sured by LSND jloj and KARMEN [ll|, agree quite well 
with each other and also with theoretical calculations. To 



reconcile the LSND, KARMEN, and MiniBooNE results 
in terms of the, so-called (3+l)-f oscillations scheme or 
(a yet unknown) experimental background seems quite 
difficult [HI]. Therefore, it is obviously important to ask 
whether neutrino oscillations are the only possible expla- 
nation for the observed anomalies. 

This work has attempted to reconcile puzzling neu- 
trino oscillation results from the LSND, KARMEN and 
MiniBooNE experiments. Our discussion is based on the 
fact that signals produced by electrons or by converted 
photons in these experiments are indistinguishable. This 
hint suggests that the excess events observed by LSND 
and MiniBooNE could originate from converted photons, 
and not from electrons. As an input, we use a natural 
extension of the model developed in Ref. 13] for an ex- 
planation of the MiniBooNe anomaly observed in data 
in terms of the radiative decays of a heavy neutrino. We 
show that the LSND evidence for — > 77 e oscillations, its 
long-standing disagreement with the results from KAR- 
MEN, and the anomalous event excess observed by Mini- 
BooNE in and data could all be explained by the 
existence of a heavy neutral lepton {yh). All these ob- 
servations are found to be consistent with each other, as- 
suming that the z/^'s are produced in neutral-current 
interactions (NC) and that they decay radiatively into a 
photon and a light neutrino v. The z/^'s could be Dirac 
or Majorana type and could decay dominantly into jis if, 
e.g., there is a large enough transition magnetic moment 
between the Vh and v mass states. Discussions of other 
decay modes suggested for the explanation of the LSND 
signal can be found in Ref. [3]. 

We may consider the Uh as a very weakly interacting 
particle directly produced by the flavor eigenstate in 
neutrino-nucleus reactions. However, it is known, that 
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the neutrino weak flavor eigenstates (y e , v^, v T , ...) can 
be different from the mass eigenstates (y\, i>2, ^3, v^...), 
but they are related to them, in general, through a uni- 
tary transformation. A generalized mixing: 
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results in neutrino oscillations when the mass differences 
are small, and in neutrino decays when the mass differ- 
ences are large. Hence, it would also be natural to assume 
that the Vh, if it exists, is a component of muon neutri- 
nos which is produced in NC interactions by muonic 
mixing, as illustrated in FigflJ This assumption provides 
us with a useful framework for further discussions. An 
immediate consequence is that the i>h can also be pro- 
duced through CC interactions in leptonic and semilep- 
tonic decays of sufficiently heavy mesons and baryons 
according to the proper mixing strength, as follows from 
Eq. ([T]), and phase space and helicity factors [H, [l6| 
(see also [D})- Note that, although CC weak interac- 
tions of ordinary particles are V — A, one could assume 
that the heavy neutrinos may dominantly be produced 
by non-left-handed V, A couplings; see e.g., the discus- 
sion in Ref.[f5]. Therefore, it would be interesting and 
important to have a general analysis of the production of 
heavy neutrinos of Dirac or Majorana type, e.g. in v^NC 
interactions, for arbitrary weak couplings including the 
leptonic mixing and helicity effects. This is, however, 
beyond the scope of the present work. 





FIG. 1: Schematic illustration of the NCQE production and 
the decay of heavy neutrino. 



The rest of the paper is organized as follows. In Sec. II 
we describe the formalism for the radiative neutrino de- 
cay, specifying the difference between the Dirac and Ma- 
jorana decay modes. The results from the LSND and 
KARMEN experiments are described in Sec. III. Here we 
show how the suggested model explains those results. In 
Sec. IV we briefly describe the MiniBooNE experiment 
and give an explanation of the anomalous excess of events 
observed in f M and data. The final results from the 
combined analysis of the LSND and MiniBooNE data are 
reported in Sec.V. The discussion and review of the ex- 
perimental and some cosmological and astrophysical con- 
straints on the mixing strength |J/^h| 2 and neutrino mag- 
netic moment are presented in Sec. VI. We find that, quite 



surprisingly, the (m„ k ; parameter space favorable 

for the explanation of the LSND and MiniBooNE results 
is unconstrained by the results from the most sensitive 
experiments, e.g. searching for a Vh peak in 71^2,-^2 de- 
cays. Moreover, we show that taking into account the 
dominance of the radiative Vh decay and its short life- 
time, makes existing experimental bounds weaker, allow- 
ing them to be extended to the higher mass region. In 
Sec. VII, several proposed experiments to search for the 
Vh are described. We also show that, several tests can be 
applied to existing data. Section VIII contains conclud- 
ing remarks. 



II. RADIATIVE NEUTRINO DECAY 



Let us consider the decay of a heavy neutrino i>h of 
mass m v , and energy E„ h into a lighter neutrino v and 



a photon: 



v + 7 



(2) 



with the partial lifetime T Vh . The energy of the decay 
photon in the Vh rest frame given by 
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is in the range < £L < m Vh /2, depending on the mass 
of the which may be in the range < m v < m Vh . 
Furthermore, for simplicity we assume that the particle 
v is almost massless, and the photon energy in the rest 
frame is E® = m Vh /2. The energy of the decay photon 
in the laboratory frame depends on the Vh initial energy 
and on the center-of-mass angle between the photon 
momentum and the Vh direction of flight: 

E -> = ^T (1 + ^7 cos6) - \ [x + cos9) (4) 

Hence, the energy distribution of photons in the labora- 
tory system depends on their angular distribution in the 
rest frame, which is not generally isotropic [Til ]: 
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^e = 2 (1 + acose) 



(5) 



Here, the angle is defined as above and a is the asym- 
metry parameter. It is also possible to define 6 as the 
angle between the direction of spin, the only direction 
available in the rest frame, and the photon momentum. 
However, if we assume that the spin of Vh is (anti)parallel 
to its momentum, both definitions are equivalent. 

The decay of a spin-^ neutrino into another spin-i 
particle and a photon can be generally described by two 
helicity amplitudes A and B corresponding to the final 
states shown in Fig. [3] For the most general coupling 
given by |19l42f| 
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FIG. 2: Two amplitudes, A and B, describing the decay 
Vh 7^ f° r different Uh helicities. 



the amplitudes A and B are proportional, respectively 
to (a — j3) and (a + f3). If CP is conserved, the helicity 
amplitudes |A| and \B\ for the decay of Majorana neu- 
trinos are equal. In this case the decay Vh — > would 
be isotropic and independent of the Vh polarization, and 
hence a = in Eq.JSJ). Indeed, suppose the vh — > jv de- 
cay is anisotropic in the center-of-mass system and pho- 
tons are emitted preferably, say opposite to the Vh spin 
direction. Because of CP conservation, the CP-mirror 
image of this process should also exist, and it would cor- 
respond to the Vh — > ji> decay with photons emitted 
preferably along the Vh spin direction. But, if the is 
its own antiparticle, the decay anisotropy must be the 
same for the Vh and its CP-mirror image. Hence, the de- 
cay must be isotropic. For the Dirac case, the Vh and its 
CP-mirror image are not identical and the above argu- 
ments do not hold. For Dirac Vh the angular anisotropy 
is the result of parity nonconservation in the decay © 
and of nonvanishing polarization of the neutrinos. The 
decay asymmetry parameter given by 
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is, in general, not constrained, and it may be in the range 
-1 < a < +1 [H, H3|. In the standard model (3 /a = 
(ra Vh — m v )/(m Uh +m v ), so that 
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m,, 



is equal to zero only when 
Dirac neutrinos and m„ <C rn v , , one has 



For left-handed 
— 1, which 

means that the decay photons are emitted preferably 
backward (l9l - [2l| . shifting the energy spectrum in the 
laboratory frame towards lower energies. For the right- 
handed Dirac neutrinos, one has a — +1, and the photons 
are emitted preferably in the forward direction, making 



the energy spectrum harder. Hence, the energy spectrum 
and angular distribution of the decay photons are sensi- 
tive to the type of Vh- Note that if CP is conserved, the 
decay rate and the center-of-mass angular distributions 
for the Dirac case are the same for the Vh — > as for the 
Vh — > V^f decay modes with respect to the beam direc- 
tion. Furthermore, we assume that the decay Vh "/v is 
generally CP conserving (see also Sec. IV. B). 

As mentioned above, the most natural way to allow 
the radiative decay of heavy neutrino is to introduce a 
nonzero transition magnetic moment (utr) between the 
Vh and v mass states; see e.g. [H, |23j. Such coupling 
of neutrinos with photons is a generic consequence of 
the finite neutrino mass. Observations of the neutrino 
magnetic moment could allow to distinguish if neutrinos 
are of Dirac or Majorana type since the Dirac neutri- 
nos can only have flavor conserving transition magnetic 
moments while the Majorana neutrinos can only have 
a changing one. In addition, Dirac neutrinos can have 
diagonal magnetic moments while Majorana neutrinos 
cannot. The nonzero magnetic moment of the neutrino, 
although tiny, is predicted even in the standard model. 
The detailed calculations of the radiative neutrino de- 
cay rate in terms of the neutrino masses and mixings of 
Eq.([T]) were performed long ago, see e.g. (l9l - l2]| . The 
radiative decay mode could even be dominant, if the pL tr 
value is large enough; see [22|, [23| . Originally, the idea of 
a large (Dirac) magnetic moment (> 10 _11 /Lts, where [Ib 
is the Bohr magneton) of the electron neutrino has been 
suggested in order to explain the solar neutrino flux vari- 
ations [1^ . Taking into account that in many extensions 
of the standard model the value of the fj,t r is typically 
proportional to the Vh mass, the intention to make the 
radiative decay of a < 100 MeV dominant by intro- 
ducing a large transition magnetic moment (or through 
another mechanism) is not particularly exotic from a the- 
oretical viewpoint. Such types of heavy neutrinos are 
present in many interesting extensions of the standard 
model, such as GUT, superstring inspired models, left- 
right symmetric models and others, for a review; see e.g. 

The total Vh decay width can be defined as T to t — 
T{vh — > vy) + Er^, where T(uh — > vj) is the Vh — > 
decay rate, and EFj is the sum over decay modes whose 
decay rate is proportional to the square of the mixing 
\U^h\ 2 - For the Vh with a mass < 100 MeV, the domi- 
nant contribution to T,Ti comes from Vh v^e^e^ and 
v^viyi (I = e, /i, r) decays, for which the rate calculations 
can be found, e.g. in |25l - l27| . The Vh — > 7^ decay rate 
due to a transition moment [itr is given by [28[ 
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The decay rate T(yh — > v^ee) can be estimated as 



T(^ -> v^e + e 1 



10 MeV 



(9) 



(10) 



and the sum rate ET^ ~ 9 • T(vh v^e + e ). For 
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50 MeV, |CV| 2 < 1CT 2 and fi tr > IO^Vb, we found 
that the radiative decay is dominant, as its branching 
fraction Br{v h -> 71/) = r( ^^ 7) > 0.99. 

III. INTERPRETATION OF THE LSND AND 
KARMEN RESULTS 

The LSND and KARMEN experiments used neutrinos 
produced in the beam stop of a proton accelerator. LSND 
finished data taking at LANSCE at the end of 1998, while 
KARMEN finished data taking in 2001. In these experi- 
ments, neutrinos were produced by the following decays 
of pions and muons occurring in the proton target: 

• 7r + —> n + v^ decays in flight (DIF) or at rest (DAR), 

• i± + ->• e + v e V^ DAR, 

• 7T~ — > fJ>~Vfj, DIF, 

• /i~ — > e~v e v^ DAR. 

The main detector properties and the neutrino fluxes in 
these experiments are summarized in Table 1. 

A. The LSND signal of z7 M — > T> e oscillations 

In 1996 the LSND experiment published evidence for 
— > ~D e oscillations, based on the observation of an ex- 
cess of F e -like events [2j ■ Measurements performed from 
1996-1998 with a different target configuration confirmed 
the evidence and improved the significance of the ob- 
served excess. The LSND detector is described in detail 
in Ref. [29[ . It was located at a distance of 30 m down- 
stream of the main LANSCE beam-stop A6 at a small 
angle of ~ 12° relative to the primary proton beam. The 
detector was a cylindrical volume filled with 167 t of a 
dilute mineral oil (CH2) based liquid scintillator viewed 
by photomultipliers (PMT) and surrounded by an ac- 
tive 47r veto shield. The low light-yield of the scintillator 
allowed for the detection of Cherenkov light generated 
by relativistic muons, electrons, and converted photon 
tracks. This feature was of great importance for parti- 
cle identification and reconstruction of its direction. The 
energy resolution of the detector was about ~ 6% at 50 
MeV electron energy. 

The search for Fp — > V e oscillations was based on the 
appearance of v e in the neutrino beam, detected through 
the reaction V e p — > e + n resulting in a prompt rela- 
tivistic e + , followed by a 2.2 MeV gamma signal from 
the neutron capture p(n,^/)d. The e + candidate events 
identification and separation from the background were 
based on the detection of the prompt and directional 
Cherenkov light, and scintillation light which is delayed 
and isotropic. The 2.2 MeV signal from the reaction 
p(n,~/)d is correlated in time with the positron one. It 
was identified and separated from accidental low-energy 



7's by means of a likelihood parameter i? 7 , which is de- 
fined as the ratio of the likelihood of a low-energy event 
being correlated or being accidental. The parameter IL, 
was defined by three values: i) the PMT multiplicity, 
which is proportional to the 7 energy, ii) the radial dis- 
tance between the reconstructed positions of the e + and 
7, and iii) the time difference between the e + and 7, 
which is defined by the capture time of 186 /is of neu- 
trons in mineral oil, while accidentals are distributed 
uniformly in time. A \ 2 fit to the i? 7 distribution ob- 
tained from the 1993-1998 measurements resulted, after 
subtraction of background from DAR and DIF neutrino 
events, (19.5 ± 3.9) and (10.5 ± 4.6), respectively, in a 
beam on-off excess of (87.9 ± 22.4 ± 6.0) events. The 
neutrino background was carefully evaluated both from 
independent measurements and calculations. This excess 
was attributed to the appearance of V e from — ¥ V e os- 
cillations and corresponds to the oscillation probability 
P{V^ -> V e ) = (2.64 ± 0.67 ± 0.45) x 10~ 3 . 



p, 800 MeV T 




i 



KARMEN 

FIG. 3: Schematic illustration of the location and orientation 
of the LSND and KARMEN detectors relative to the incident 
proton beam direction. 



The KARMEN experiment used a technique similar to 
the LSND experiment and observed no beam excess 
The signatures of the 15 candidate events were found to 
be in good agreement with those from the (15.8 ± 0.5) 
expected background events. 

Let us explain the discrepancy between the results of 
these two experiments in terms of the production and ra- 
diative decay of a heavy neutrino, as illustrated in FigJT] 
The location of the LSND and KARMEN detectors rel- 
ative to their proton beam directions is schematically il- 
lustrated in Fig. |31 The energy distributions of the v^s 
from 7T+ DIF in the LSND (from Ref.[H) and KARMEN 
(simulated) detectors are shown in FigfJJ The distribu- 
tions are normalized to a common maximum value in 
order to place them on a similar scale. One can see that 
the spectra are quite different. The LSND distribution is 
peaked at about 55 MeV; it has an average energy ~ 100 
MeV, and a high energy tail up to ~ 300 MeV, while the 
maximum of the energy spectrum in KARMEN, which is 
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TABLE I: Comparison of experimental parameters of the LSND and KARMEN experiments. 



LSND 



KARMEN 



p beam kinetic energy, MeV 
total number of POT's 
distance to target, m 
angle between the v and p beams 
total 77^ flux 
77 M ,^ e /cm 2 from fi + DAR 
^ M ,77 e /cm 2 from fi~ DAR 
v^/cm 2 from tv + DIF 
VcP —5- e + n efficiency 
e + energy range , MeV 
observed events 
background 
event excess, i? 7 > 1 
event excess, J? 7 > 10 
77 M — ► V e oscillation probability 



800 
1.8 x 10 23 
30 
12° 
1.2 x 10 22 
1.26xl0 14 
1. 08x10" 
2.2xl0 12 
0.17 
20 - 60 
87 
53.8 

87.9 ±22.4 ±6.0 
32.2 ±9.4 
(2.64 ±0.67 ±0.45) x 10~ 3 



800 
5.9 x 10 22 
17 
90° 
2.71 x 10 21 
8.86xl0 13 
7.6xl0 10 
< 10 11 

0.19 
16 - 50 
15 

15.8 ±0.5 
10 ± 32 
< 5.1(90% C.L.) 
< 0.85 x 10" 3 (90% C.L. 



located at 90° with respect to the beam, is ~ 20 MeV and 
the whole spectrum is well below 50 MeV. For a heavy 
neutrino with a mass of m Vh = 40 MeV, the production 
threshold in the reaction vjfC — > v h n 11 C g . s . is 58.6 MeV, 
as shown in FigJU Here we assume that the Vh produc- 
tion is accompanied by the emission of a recoil neutron 
and the isotope 11 C in the ground state. 

Thus, our interpretation of the excess of events ob- 
served by LSND is the following. Positive pions gener- 
ated in proton collisions produce the flux of z^'s from 
the 7r + — » /Lt + i^ t DIF in the target. The excess events are 
generated in the LSND detector by these v^s through 
the reaction 

vfC -> v h nX -> "fvnX, (11) 

with the emission of a recoil neutron and a heavy neu- 
trino, and not by 77^'s from muon decays at rest via 
Vfj — > v e oscillations, as was originally assumed Q. The 
VhS, decay promptly into a photon and a light neutrino, 
with the subsequent Compton scattering or e + e _ pah- 
conversion of the decay photon in the detector fiducial 
volume. The former process dominates for photon ener- 
gies below the critical energy of the LSND liquid of 85 
MeV. In the laboratory system, the differential Compton 
scattering cross section has a sharp peak in the forward 
direction, and the vast majority of events are in a nar- 
row cone of < y Jm e /E 1 < 100 mrad, for E 1 > 20 MeV. 
For the photon conversion into an e + e~ pair, its opening 
angle is ~ m e /E 1 < 25 mrad for E 1 > 20 MeV, which 
is too small to be resolved in LSND into two separate 
Cherenkov rings (here, m e ,iJ 7 are the electron mass and 
the photon energy, respectively). Therefore, the excess 
events are originated from photons of the reaction (| 1 1|) 
detected in coincidence with the associated 2.2 MeV 7 
tag from the neutron capture and misidcntified as single 
electron events. In the KARMEN experiment, u^s from 
7T decays in flight cannot produce heavy neutrinos accom- 
panied by the emission of a neutron because their energy 
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FIG. 4: The shape of the energy distributions of the u^s 
from 7T+ DIF in the LSND (9 = 12°) and KARMEN (9 = 
90°, hatched) detectors. The arrow shows the production 
threshold of E t h = 58.6 MeV for the heavy neutrino with a 
mass of 40 MeV in the reaction v ^ ± 12 C Vh + n + 11 C g . a ., 
in which the Vh production is accompanied by the emission 
of a neutron and the isotope 11 C in the ground state. The 
distributions are normalized to a common maximum value. 



is below the Vh production threshold, see Fig. 2) There- 
fore, KARMEN should observe no excess of 77^ — > v e - 
like events. Note that the maximum energy of F^'s from 
muon DAR is about 50 MeV and is also less than the 
energy threshold of 58.6 MeV for the production of the 
40 MeV Vh and a recoil neutron in collisions with the 
carbon nucleus. 

To make quantitative estimates, we performed simpli- 
fied simulations of the Vh production in the inclusive re- 
action (jlip with the emission of a recoil neutron and fol- 
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lowed by the decay — > jv, as shown in Fig. [T] in the 
LSND detector. In these simulations we used the integral 
DIF energy spectrum, shown in Fig. |4j which was cal- 
culated in [30] • There is also a contribution from F M DIF 
events, which, however, is small and is neglected at the 
level of accuracy of our analysis. The energy of most of 
the v^'s is well above the threshold for the production of 
40-80 MeV v h 's in the LSND detector. Once produced, 
the fh's decay at an average distance ~ cThE Uh /m Uh from 
the primary vertex. Since in the LSND experiment the 
average Vh kinetic energy is E Uh ~ 50 MeV and zVs 
would decay over the average distance of < 5 m from 
the primary vertex, the sensitivity is restricted to the Vh 
lifetimes T Uh < 10~ 8 s for the Vh masses m Vh > 40 MeV. 
The decay photon absorption occurs at a distance of the 
order of the Compton scattering length (~ 40 cm ) of the 
LSND liquid from the Vh decay point which is much less 
than the detector size. 




10 20 30 40 50 60 70 80 90 100 
Neutron kinetic energy, MeV 



FIG. 5: The shape of the distributions of the recoil neu- 
tron kinetic energy in the reaction f„ C — > v^nX obtained 
in the present work (solid histogram ) and calculated in [3l|] 
(solid curve ) for the v^ energy E v =150 MeV. The distribu- 
tion of the kinetic energy of neutrons ejected in the reaction 
v^C — > VhnX for a heavy neutrino mass of 60 MeV is also 
shown for comparison (dashed histogram). The binding en- 
ergy corrections are not applied. The distributions are nor- 
malized to a common maximum value. 



The total cross section of the reaction vj^C — > vyjiX 
for 100% mixing is estimated by extrapolating the avail- 
able cross section for the reaction vj^C — > v^nX (~ 
24 x 10 -40 cm 2 ) calculated for the incident neutrino en- 
ergy of 150 MeV [111, HI] to the neutrino energies in the 
range 50-250 MeV (see Fig. 0}, an d by taking into ac- 
count the corresponding phase space factor. Note that 
the average energy for the spectrum above the pro- 
duction threshold of the 40 MeV v h is 110 MeV. The 
Vy, flux-averaged cross section is found to be u(v 12 C — > 
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2 3000 - 




10 20 30 40 50 60 70 80 90 100 
Neutron kinetic energy, MeV 

FIG. 6: The flux-averaged distributions of the recoil neu- 
tron kinetic energy from the reaction v^ 2 C — >• v^nX calcu- 
lated for Vh masses of 40 and 80 MeV with the binding en- 
ergy correction included. The histograms are normalized as 
in FigOU 



VhnX) ~ (16±6.5) x 10 -40 cm 2 for the production of the 
40 MeV Vh with the uncertainty taken to be 40% due to 
accuracy of the extrapolation procedure. 

The crucial test of the v^ — > v e oscillation hypothesis 
in the LSND experiment was to check whether there is an 
excess of events with more than one correlated 2.2 MeV 
7. If the excess of events is indeed due to the reaction 
v e p — > e + n, then there should be no excess with more 
than one correlated 7 because the recoil neutron is too 
low in energy ( < 5 MeV) to knock out additional neu- 
trons. If, on the other hand, the excess involves higher 
energy neutrons, which can break the 12 C nucleus and 
produce another neutron(s), then one would expect an 
excess of events with > 1 correlated 7. As the LSND 
did not observe many of such (latter) events the en- 
ergy spectrum of the ejected neutrons is an important 
characteristic of the reaction (ITTT) . as it affects the likeli- 
hood ratio i? 7 and the number of correlated 7's. The v^ 
flux-averaged cross sections and particle emission spec- 
tra in the LSND detector are predicted quit well for the 
charge current reactions v^ 2 C — > uX from calculations 
performed in different theoretical frameworks; see e.g. 
Ref.[33j and references therein. However, much less is 
known for the v^ induced neutral-current reactions in the 
detector. Therefore, we performed simulations of the re- 
coil neutron kinetic energy distributions of in the reaction 
(|11[) by using the Fermi gas nuclei model, but without 
taking into account nuclear effects, such as the neutron 
re-scattering in nuclear matter and the carbon nucleus 
level structure. The Fermi momentum and the neutron 
binding energy for the 12 C nucleus are taken to be 200 
MeV and 18 MeV, respectively. 



7 



> 
CD 



10000 



8000 
° 6000 
4000 
2000 




_Q 




10000 
8000 
6000 
4000 
2000 





50 100 150 



50 100 150 



> 6000 

CD 

^ 5000 

}r 4000 

2 3000 

CD 

f 2000 
z 1000 




,a = -1 



m„ h = 70 MeV 



a = 




50 100 150 
E r MeV 




50 100 150 
E r MeV 



FIG. 7: The flux-averaged distributions of the energy of photons from the radiative decay of heavy neutrinos produced in 
the reaction v x ^C — > v^nX calculated for 100% mixing strength and v% masses from 40 to 80 MeV with no photon detection 
efficiency and the neutron binding energy correction included . The spectra are calculated for the a = — 1 (solid line), a = +1 
(dotted line), and a — (dashed line) cases for the same i/ M flux. 



To evaluate uncertainties of our calculations we have 
compared our results with others which take into account 
nuclear effects 31]-35|. Figure^] shows the distribution 



of the kinetic energy of neutrons ejected in the reaction 
v^C — > v^nX and the analogous spectrum from Ref . [3l| . 
both calculated for the massless case for the incident neu- 
trino energy of 150 MeV without nucleon binding energy 
corrections. One can see that our simulations reproduce 



the more precise results quite reasonably. The compari- 
son of the calculations results in an uncertainty of about 
20%-30%. Figure|S] also shows the neutron energy distri- 
bution calculated for the reaction (fTT|) for the Vh mass of 
60 MeV. It is seen that in this case the neutron energy 
spectrum is shifted towards lower energies. 

Having this reasonable agreement in mind, we have 
performed calculations of the LSND flux-averaged dis- 
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FIG. 8: Distributions of the excess events reconstructed as z7 e CC events in the LSND detector as a function of visible energy 
E V i S for R-y > 10 from the 1993-1998 data sample (dots), and from a combination of the — > decay plus expected neutrino 
background calculated for a = — 1 (solid line) and a — +1 (dashed line), vu masses of 40 and 70 MeV shown in the plots, the 
mixing strength |£/ M h| 2 = 3 x 10 -3 , and the vu lifetime r Vh = 10 -9 s. A combination of neutrino background plus neutrino 
oscillations at low Am 2 (dotted line) and the expected distribution from neutrino background (shaded, from Ref. [3]) are also 
shown. 
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FIG. 9: The same as Fig. [8] for the case a = 0. 



tributions of the kinetic energy of knockout neutrons pro- 
duced in the reaction (fTTj) . The results are shown in Fig IB] 
for the Vh masses of 40 and 80 MeV. The average energies 
of the recoil neutrons are, respectively, 14 and 16 MeV. 
To decrease this energy to the typical energy of neutrons 
from the reaction T> e p — > e + n (< 5 MeV) takes about 
n co u ~ 6 neutron collisions in mineral oil. Taking into 
account that the average collision length is L co u < 10 cm 
results in a displacement of the neutron from the primary 
vertex of the order Ar ~ L co u^/n co u ~ 25 cm, which is 
significantly less than the value of ~ 70 cm for the neu- 
trons from the reaction v e p — ¥ e + n, defined mainly by the 



reconstruction accuracy [3]. Thus, one would expect no 
significant contribution to the likelihood ratio i? 7 due to 
this effect. The energy decreasing time At ~ Ar/ (3c ~ 5 
ns is also much less than the neutron capture time of 186 
/j,s. 

As discussed above, the neutrons from the higher en- 
ergy tail of the distribution shown in Fig. [5] can knock 
out an additional neutron(s), resulting in the observa- 
tion in LSND of a number of events with more than one 
correlated capture 7's. To estimate this number, we use 
the results of the measurement of the neutron yield from 
the 70 MeV proton beam collisions with a thick graphite 
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target [36], assuming that this yield is approximatelly 
the same for the neutron induced reaction of the same 
energy. The measured number of neutrons per proton 
is found to be ~ 0.06. The neutron energy threshold 
to produce a secondary neutron in collisions with a 12 C 
nucleus is about 20 MeV. The fraction of neutrons with 
energy greater than 20 MeV in the distribution shown in 
Fig|51 is ~ 20 — 25% depending on the value of the Vh 
mass. Taking this into account, we find that the total 
fraction of events with more than one correlated gamma 
from the flux-averaged reaction ([TT|) is < 2%, which can 
be neglected. 

Note that in our calculations we overestimate the frac- 
tion of high energy neutrons. The calculations of the 



reaction v¥C — > fi v n 11 C performed in Ref.(34| show 
that the cross section and recoil neutron energy spec- 
trum are essentially dependent on the details of the 12 C 
level structure for neutron energies below 30 MeV. The 
consideration of such an effect, including the rescattering 
of outgoing neutrons, is quite important for the emission 
spectrum, as it will shift the spectrum to lower neutron 
energies. Therefore, we may assume, that the fraction 
of events with > 1 correlated 7 is even less than 2%, or 
0.6 events. This number should be compared with the 
background of < 5 such events expected in the LSND ex- 
periment at the 2 a level ||. Thus, our estimate is com- 
patible with the (approximately zero) number of events 
with > 1 correlated 7's observed by LSND for the full 
20 < E ms < 60 MeV energy region. 

A cross-check for the fraction of neutron emitted in the 
reaction (|lip can be obtained from the comparison of re- 
sults obtained for the inclusive reaction 12 C(F M , /j, + n)X 
in this work, in Ref.(33[, and by the LSND Collabora- 
tion |30| | . In this reaction the presence of a muon and 
a neutron is established by detection of the Chercnkov 
ring and of the 7 ray from the neutron's capture as de- 
scribed in detail in Ref. [13] ■ We found that the fraction 
of events accompanied by the emission of a recoil neu- 
tron is ~ 81%. This number has to be compared with 
the fractions of 79 % predicted from the calculations in 
33], and (79.6 ± 40.0)% obtained by LSND [30|. The 
agreement is quite good. 

In Fig [7] the energy distributions of photons from the 
radiative decay of heavy neutrinos produced in the re- 
action v 12 C — > v^nX calculated for several Vh masses 
from 40 to 80 MeV and for three of the most interest- 
ing cases of the decay asymmetry parameter, a = ±1 
and a = 0, are shown. The calculations are performed 
with no photon efficiency and no binding energy correc- 
tions included. One can see that for the Dirac case with 
a = — 1, the simulated events are mainly distributed in 
the narrow region < E 1 < 60 MeV. The fraction of 
photons in this region varies from 0.86 to 0.77 for a i>h 
mass from 40 to 80 MeV, respectively. The remaining 
events are distributed over the region 60 < E 1 < 150 
MeV, where they can be hidden by the low statistics. 
For the a = and a = +1 cases, the fraction of photon 
events in the region < E 1 < 60 MeV varies from 0.71 



to 0.57; and from 0.54 to 0.46, respectively. Further, we 
will discuss mainly the cases a = — 1 and a = 0, because 
for the case a = +1 the fraction of events above 60 MeV 
is too high compared to the LSND observations. 

The visible energy distributions expected from a com- 
bination of the Vh — > 71/ events plus neutrino background 
in the LSND detector Q are shown in Figs. 18191 for the 
energy range 20 < E V i s < 60 MeV. The spectra are cal- 



culated for 



40 and 70 MeV, \U, 



3x10-3, and 



the Vh lifetime r Vh — 10 s by taking into account the 
decay photon efficiency and the neutron binding energy 
corrections. The photon efficiency has been estimated by 
a simple Monte Carlo calculation as a fraction of photons 
with energies 20 < E y < 60 MeV in the detector fiducial 
volume times the detection efficiency, which is taken to 
be essentially constant, e 7 ~ 0.4, for this energy range. 
The distribution expected from a combination of neu- 
trino background plus neutrino oscillations at low Am 2 
Q is also shown for comparison. The clean experimental 
sample of the oscillation candidate events shown in Fig.'s 
is obtained by enforcing strongly correlated gammas 



with the cut i? 7 > 10. In this case the beam on-off 
excess is 49.1 ± 9.4 events while the estimated neutrino 
background is only 16.9 ± 2.3 resulting in a total excess 
of 32.2 ±9.4 ±2.3 events [|[. The analogous distributions 
for cos0 7 „, the cosine of the angle between the incident 
neutrino beam and decay photon momenta, for events 
with 36 < E ms < 60 MeV are shown in Figs. US HU 
The distributions are obtained assuming that the energy 
deposited by the decay photon is misreconstructed as the 
energy from a single electron track. 

Simulations are in reasonable agreement with the ex- 
perimental distributions. For instance, for the distribu- 
tion shown in Fig|S] for m Uh = 40 MeV, the comparison 
with the LSND data yields a x 2 of 3.6 for 3 DF, corre- 
sponding to 34% C.L. The best fit results suggest that 
the Vh mass is in the region 10 < m Uh < 90 MeV and the 
lifetime is t v% < 10~ 8 s. However, to avoid the produc- 
tion of Vh's in the KARMEN experiment, the low mass 
limit is set to 40 MeV. The mass upper bound is set to 
80 MeV because for higher masses the production of Vh 
in LSND is suppressed by the phase space factor. The 
simulations showed that the shape of the E V i S distribu- 
tion is sensitive to the choice of the Vh mass : the higher 
the mass, the harder the visible energy spectrum. 

Before the calculation of the required mixing strength 
Wfih] 2 : let us estimate, for a cross-check, the number 
of events expected for v^ — > v e oscillations followed by 
v e p — > e + n scattering in the LSND detector. This num- 
ber could be estimated as 



A§ v P osc <j v J e z e 



(12) 



where A — 7.4 x 10 30 is the number of free protons in the 
LSND fiducial volume, $77 is the neutrino flux 1.26 x 
10 14 v/cm 2 (see Table 1), P osc is the v^ — > v e oscillation 
probability averaged over the incident neutrino energy, 
0> e = 0.95 x 10~ 40 cm 2 is the cross section averaged 
over the entire energy range, f e ~ 0.89 is the fraction 
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FIG. 10: Distributions of the excess events reconstructed as V e CC events with 36 < E V i„ < 60 MeV in the LSND detector as 
a function of cosO from the 1993-1998 data sample (dots), and from a combination of the the Vh jv decay plus expected 
neutrino background calculated for a = — 1 (solid line) and a = +1 (dashed line), vu masses of 40 and 70 MeV shown in the 
plots, the mixing strength |t/ M h| 2 = 3 x 10 -3 , and the Vh lifetime r Uh = 10 -9 s. A combination of neutrino background plus 
neutrino oscillations at low Am 2 (dotted line) and the expected distribution from neutrino background (shaded, from Ref. [3]) 
are also shown. 
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FIG. 11: The same as in Fig. QJJfor the case a = 0. 



of events in the energy range 20 < E < 60 MeV, and 
e = 0.42 is the average positron reconstruction efficiency 
Q. Using the above values, we found that the LSND 
experiment should detect an excess of AA^ _>i7 e — 70 
events if the oscillation probability is P osc ~ 2.6 x 10~ 3 . 
This value is in a good agreement with the number of 
events 87.9 ± 22.4 ± 6.0 quoted in Ref.[l]. 

Consider now the case of heavy neutrino. The estimate 
of the mixing parameter lET^I 2 was performed by using 
the following relations. For a given flux of muon neutri- 
nos, <&„ , the expected number of the Vh —> 7^ decays 



events in the LSND detector is given by 

AN Vh ^ v ~ A J $„ M Vfk | U^hfnhhsPdeoPabs^dE 

(13) 

where £±N Vh ^ lv — 32.2 ± 9.7 is the number of excess 
events observed in the 1993-1998 data sample (with er- 
rors combined in quadrature), A = 3.7 x 10 30 is the num- 
ber of carbon nuclear in the LSND fiducial volume, <J Viin 
is the cross-section for the reaction 12 C(f /i , v^n)X with 
the emission of a recoil neutron for the massless case, 
/ 7 ~ 0.5 — 0.76 is the fraction of events in the energy 
range 20 < E < 60 MeV, /„ ~ 0.4 - 0.8 is the fraction 
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of events with the emission of a recoil neutron in reac- 
tion (jlll) . Pfiec is the probability for the — > decay 
within the detector fiducial volume, P a b s is the probabil- 
ity of the decay photon absorption in the detector, and 
e 7 ~ 0.4 is the overall efficiency for decay photon detec- 
tion. In Eq.(|13p the number of heavy neutrinos produced 
is proportional to the product of the 12 C(^ M , v^n)X cross 
section, the mixing |C/ M h| 2 , and the phase space factor 
fphs, which takes into account the threshold effect due to 
the heavy neutrino mass. The value of the total f„ DIF 
flux (see Table 1), the number of reconstructed — s- V e - 
like [3[ and v^CC events [30| in the detector were used 
for cross checks and normalization. The flux averaged 
12 C(u^, v^njX cross section is < a(E)y^ n >~ 16 x 10~ 40 
cm 2 . The probability of the heavy neutrino to decay ra- 
diatively in the fiducial volume at a distance r from the 
primary vertex is given by 

^ c = [l-ex P (^)]£^^ (14) 

Pu h T Uh I tot 

where the last term is the branching fraction Br{yh — > 
71/) ~ 1. Assuming that almost all — ¥ 71/ decays occur 
inside the fiducial volume of the detector, we estimate the 
\Ufi.h\ 2 to be in the range 

|t^fc| 2 ~ (3-9) x 1(T 3 . (15) 

This result is mainly defined by the uncertainty on the 
number of excess events and is valid for the mass region 

40 < m Vh < 80 MeV (16) 

and the lifetime 

Tu h < 11T 8 « (17) 

B. The LSND signal of — s> v e oscillations 

During the first three years of LSND data taking, the 
target area of the LANSCE accelerator consisted of a 
30 cm long water target located ~1 m upstream of the 
beam stop. This configuration enhanced the probabil- 
ity of pion decay in flight, allowing LSND to search for 

— v e oscillations using with energy above 60 MeV. 
In this case, one expects to observe an excess of events 
from the reaction v\ 2 C — > e~X above the expected back- 
grounds. This reaction has only one signature (a prompt 
signal), but the higher energy and the longer track of 
the events allow good electron identification and measur- 
ing its direction. In this search LSND has observed 40 
events to be compared with 12.3 ± 0.9 events from cos- 
mic ray background and 9.6 ± 1.9 events from machine- 
related (neutrino-induced) processes 38]. The excess of 
(18.1 ± 6.6) events corresponds to a — v e oscillation 
probability of (2.6 ± 1.0) x 10~ 3 , consistent with the 
value found from the study of the V e p e + n reaction 
below 60 MeV. 



The number of the Vh — ► 7^ events that would be 
observed by LSND after applying the high energy cut 
E > 60 MeV, is about 3 — 10 events depending on the vu 
mass and mixing obtained from the combined analysis, as 
shown below in Sec.V. For example, out of 10 events, ~ 5 
(~ 2) events are from the reaction ic 2 C — > ^hX occurring 
on protons (neutrons) of the 12 C nucleus, which is not 
expected to produce free neutrons, and ~ 3 events are 
from the reaction ((TTj) with a recoil neutron production, 
which is identified by the presence of the 2.2 MeV photon 
from the capture reaction. Thus, the ratio of the number 
of excess events with and without photon tag is ~ 3 : 
7, which is in agreement within errors with the observed 
numbers ~ (4 ± 2.5) : (15 ± 5) of events in the LSND 
experiment [38j . 

IV. THE MINIBOONE ANOMALIES 

The MiniBooNE detector is described in detail in Ref. 
[39l | . It uses an almost pure beam originated from 
the 7r + decays in flight, which are generated by 8 GeV 
protons from the FNAL booster. The detector consists 
of a target, which is a 12.2 m diameter sphere filled with 
800 t of mineral oil, surrounded by an outer veto region. 
The Cherenkov light rings generated by muon, electron 
and converted photon tracks are used for the reconstruc- 
tion of the events. The resolutions reached on the vertex 
position, the outgoing particle direction and the visible 
energy are 20 cm, 4°, and 12%, respectively for CCQE 
electrons [iol ]. The beam is peaked around ~ 600 
MeV, has a mean energy of ~ 800 MeV and a high en- 
ergy tail up to - 3 GeV [Hj]. 

Below, we consider the MiniBooNE anomalous event 
excess observed in and data and the interpretation 
of these results in terms of the heavy neutrino decay. 

A. Excess of events in data 

An excess of AJV =128.8±20.4 ± 38.3 electronlike 
events has been observed in the data accumulated with 
6.64 x 10 20 protons on target. For the following dis- 
cussion several distinctive features of the excess events 
are of importance a) the excess is observed for sin- 
gle track events, originating either from an electron, or 
from a photon converted into a e + e~ pair with a typi- 
cal opening angle ~ m e /E e + e - < 1° (for E e+e - > 200 
MeV) , which is too small to be resolved into two separate 
Cherenkov rings (here, m e ,E e + e - are the electron mass 
and the e + e~ pair energy); b) the reconstructed neutrino 
energy is in the range 200 < E$ E < 475 MeV, while there 
is no significant excess for the region E® E > 475 MeV 
(the variable E® E is calculated under the assumption 
that the observed electron track originates from v e QE 
interaction); c) the visible energy E vis is in the narrow 
region 200 < E vis < 400 MeV for events with E$ E > 200 
MeV; and d) the angular distribution of the excess events 
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FIG. 12: Distributions of the excess events in the MiniBooNE 
detector from the — » jv decay reconstructed as v e CC 
events as a function of E® E for |(7 M h| 2 = 3 x 1CT 3 , m„ h = 40 
MeV, and r„ h = 10 -9 s, and for different values of the asym- 
metry parameter a. The dots are experimental points for the 
excess events in the MiniBooNE detector. Error bars include 
both statistical and systematic errors 0. The comparison of 
the distributions with the experimental data yields a \ 2 of 
7.1 (a = -1), 9.3 (a = 0), and 10.1 (a = +1) for 8 DF. 



with respect to the incident neutrino direction is wide 
and consistent with the shape expected from v e CC in- 
teractions. To satisfy the criteria a)-d), we propose 
that the excess events are originated from the decay of 
the heavy neutrino Vh considered in Sec. III. The ish's are 
produced by mixing in neutral-current QE interac- 
tions and depositing their energy via the visible decay 
mode Vh vj, as shown in FigJTJ with the subsequent 
conversion of the decay photon into the e + e~ pair in the 
MiniBooNE target. To make a quantitative estimate, we 
performed simplified simulations of the production and 
decay processes shown in FigjT] In these simulations we 
used a energy spectrum parametrized from the recon- 
structed v^CCQE events [|l|. Since in the MiniBooNE 
experiment the z/^'s have higher energies and decay over 
an average distance of < 5 m from the production vertex, 
the sensitivity in the LSND Vh mass range of Eq. (fl~6)) is 
restricted to the lifetimes 

Tu h < l(T 9 s, (18) 

to be compared with (fl~7f . 

In Figs fP2lfT71 the distributions of the kinematic vari- 
ables E® E , E V i S and cos 71/ for the Vh — > 7^ events are 
shown for m Vh — 40 and 70 MeV and T Vh = 10~ 9 s. These 
distributions were obtained assuming that the e + e~ pair 
from the converted photon is misreconstructed as a single 
track from the v e QE reaction. Simulations are in rea- 
sonable agreement with the experimental distributions. 
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FIG. 13: Same as Fig[l2]for the 70 MeV v h . The comparison 
of the distributions with the experimental data yields a x 2 of 
7.5 (a = -1), 9.2 (a = 0), and 10.3 (a = +1) for 8 DF. 



For instance, for the distributions shown in Figs ll2U"3l 
for the case a = — 1, the comparison with MiniBooNE 
data yields a x 2 of 7.1 (7.5) for 8 DF corresponding 
to 53% (~ 47%) C.L. for m Vh = 40 (70) MeV and 
T Vh = 10~ 9 s. The simulated excess events, shown in 
Figs ll4lT5l are mainly distributed in the narrow region 
200 < E vis < 400 MeV. Their fraction in this region is 
~ 70%. The remaining events are distributed over the re- 
gion 400 < E ms < 1200 MeV, where they can be hidden 
by the low statistics. 

The simulations show that the shape of the E^ E and 
E V i S distributions is sensitive to the choice of the Vh mass: 
the heavier the Vh, the harder the visible energy spec- 
trum. The best fit results suggest that the Vh mass is in 
the region 20 < m Vh < 600 MeV and the lifetime is in 
the range r Uh < 10~ 9 — 10~ 7 s, respectively; see also []~3j j . 

The estimate of the mixing parameter |[/ M h| 2 was per- 
formed by using a relation similar to Eq. (|13[) . The flux 
$(i//j) was estimated from the expected number of the 
v^NC events times the mixing IL^I 2 , taking into ac- 
count the threshold effect due to the heavy neutrino 
mass. The total number of reconstructed v^CC events 
in the detector [4l[ was used for normalization. The 
probability of the heavy neutrino to decay radiatively 
in the fiducial volume at a distance r from the pri- 
mary vertex is given by Eq. (|14l) . assuming the branch- 
ing fraction Br(vh 7&0 — 1- Taking into account 
the ratio v^NCQE /v^CCQE ~ 0.43, and the number 
of v^CCQE events observed @, H] and assuming that al- 
most all Vh ~* l v decays occur inside the fiducial volume 
of the detector, we estimate the It^l 2 to be in the range 

|(7^| 2 ~(l-4)xl0^ 3 . (19) 

This result is mainly defined by the uncertainty on the 
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FIG. 14: Distributions of the excess events in the Mini- 
BooNE detector from the Vh — ¥ ju decay reconstructed as 
v e CC events as a function of E vis for E® E > 200 MeV, 
\U„ h \ 2 = 3 x 10" 3 , m Vh = 40 MeV and r Vh = 10" 9 s, and for 
different values of the asymmetry parameter a. The dots are 
experimental points for the excess events in the MiniBooNE 
detector. Error bars include both statistical and systematic 
errors Q]. The comparison of the distributions with the ex- 
perimental data yields a \ 2 °f 7.2 (a — —1), 9.4 (a — 0), and 
10.3 (a = +1) for 8 DF. 



number of excess events. Equation ([19]) is valid for the 
mass region 40 < m Uh < 80 MeV favored by the LSND 
data. 



B. Excess of events in data 

Recently, the MiniBooNE experiment has reported re- 
sults from a search for v ^ — \ V e oscillations using a data 
sample corresponding to 5.66 x 10 20 protons on target 
9]. An excess of AN =43.2±22.5 electronlikc events 
is observed which, when constrained by the observed 
events, has a probability for consistency with the 
background-only hypothesis of 0.5% in the oscillation- 
sensitive energy range of 475 < E < 1250 MeV. The 
data are consistent with — > V e oscillations in the 0.1 
eV range and with the evidence for antineutrino oscilla- 
tions from the LSND. Note, that the low statistics an- 
tineutrino data collected by the MiniBooNe experiment 
seem to show no low-energy excess [42[ . 

Similar to the neutrino data 0: a) the excess is ob- 
served for single track events, originating either from an 
electron, or from a photon converted into an e + e~ pair 
; b) the reconstructed neutrino energy is in the wider 
range 200 < E® E < 700 MeV, and there is also an ex- 
cess for the region E® E > 475 MeV. The variable E® E is 
calculated under the assumption that the observed elec- 



FIG. 15: Same as FiglTHfor the 70 MeV v h . The comparison 
of the distributions with the experimental data yields a x 2 of 
6.8 (a = -1), 8.8 (a = 0), and 9.7 (a = +1) for 8 DF. 



tron track originates from a u e interaction; c) compare to 
the Up data the visible energy E V i S is in the wider range 
200 < E ms < 700 MeV for events with E$ E > 200 MeV; 
and d) the angular distribution of the excess events with 
respect to the incident neutrino direction is wide and 
consistent with the shape expected from v e CC interac- 
tions. To satisfy the criteria a)-d), we propose again that 
the F M excess events are originated from the decay of the 
heavy neutrino considered in the previous sections. The 
f/j's are produced by mixing in F M NCQE interactions 
and deposit their energy via the radiative decay mode, 
as shown in FiglTJ with the subsequent conversion of the 
decay photon into an e + e~ pair in the MiniBooNE tar- 
get. 

For simulations we used a energy spectrum 
parametrized from the reconstructed v^CCQE events 
|4l| . Note that the antineutrino energy distribution has a 
maximum at ~ 400 MeV and an average energy of about 
600 MeV. 

In Figs ITBll2"3l the distributions of the kinematic vari- 

7^ events are 
= 10~ 9 s. These 



ables E® E 7 E vis and cos0 7l/ for the i>h 
shown for m Vh — 40 and 70 MeV and 7v 
distributions were obtained assuming that the e + e~ pair 
from the converted photon is misreconstructed as a sin- 
gle track from the T> e QE reaction. In this calculation we 
assume that the angular distribution of photons in the 
Vh rest frame has the same asymmetry as for the case 
due to CP conservation; see Sec. II. 

Simulations are in reasonable agreement with the ex- 
perimental distributions. For instance, for the E® E - 
distributions shown in Figs ll8ll"9l for the Dirac case with 
a = — 1, the comparison with MiniBooNE data yields a 
X 2 of 8.2 (8.3) for 9 DF corresponding to 47% (~ 45%) 
C.L. for m Vh = 40(70) MeV and r„ h = 10~ 9 s. The E ms 
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FIG. 16: Distribution of the excess events in the Mini- 
BooNE detector from the decay reconstructed as v e CC 
events as a function of cos 7 „ for 300 < E$ E < 400 MeV, 
\U^ h \ 2 = 3 x 10" 3 , m Vh = 40 MeV, and r Vfi = 10" 9 s, and for 
different values of the asymmetry parameter a. The dots are 
experimental points for the excess events in the MiniBooNE 
detector. Error bars include both statistical and systematic 
errors The comparison of the distributions with the ex- 
perimental data yields a \ 2 of 10.1 (a — —1), 9.8 (a = 0), and 
9.7 (a = +1) for 8 DF. 



FIG. 17: Same as Fig[l6]for the 70 MeV u h . The comparison 
of the distributions with the experimental data yields a \ 2 of 
9.3 (a = -1), 10.1 (a = 0), and 10.0 (a = +1) for 8 DF. 



which is consistent with the mixing from Eq. (|15p . This 
result is mainly defined by the uncertainty in the number 
of excess events. Equation (j2"0"f is valid for the mass region 
40 < m Vh < 600 MeV, which includes the region favored 
by the LSND data. 



distributions shown in Figs. 1201211 are also in reasonable 
agreement with the experiment. For the case a = — 1, the 
comparison with data yields a % 2 of 9.5 (8.5) for 5 DF cor- 
responding to 10% (~ 14%) C.L. for m Vh = 40(70) MeV 
and T Vft = 10~ 9 s. The events are mainly distributed in 
the region 200 < E V i S < 600 MeV, where their fraction 
is ~ 90%. The remaining events are distributed over the 
region 600 < E V i S < 1200 MeV, where the observed num- 
ber of events is found to be consistent with the expected 
one. 

The analysis of these data within the framework dis- 
cussed above suggests that a smaller excess of events is 
observed mainly due to the lower 77^ energy and NC cross 
section. The estimate of the mixing parameter iL^^j 2 was 
performed by using relation similar to Eq. (|13p , assuming 
the branching fraction Br(vh —¥ jv) — 1. The flux Q(vh) 
was estimated from the expected number olv^NC events 
times the mixing |{7^/,| , taking into account the phase 
space factor due to the heavy neutrino mass. The total 
number of 27,771 reconstructed v^CCQE events in the 
detector Q was used for normalization. Taking into ac- 
count the ratio V^NCQE/V^CCQE ~ 0.41 and the total 
number oiv^CCQE events observed, and assuming that 
almost all vh — > decays occur inside the fiducial vol- 
ume of the detector, we estimate the |[^h| 2 to be in the 
range 

|[/^| 2 ~(0-8)xl0- 3 . (20) 



V. RESULTS FROM THE LSND AND 
MINIBOONE DATA AND PROPERTIES OF THE 
HEAVY NEUTRINO 

To obtain the combined regions in the {m Vh ; \U^h\ 2 ) 
parameter space, we have analyzed both, the LSND 
and MiniBooNe data simultaneously. The E V i S and 
cos<d lu distributions from LSND, and the E® E , E V i S , 
and cos0 7 „ distributions from MiniBooNE were used for 
comparison with the corresponding simulated distribu- 
tions from the i>h — > decay to constrain the mixing 
strength and heavy neutrino mass. The shape of the sim- 
ulated distributions is defined by the mass ( and type) 
of the fh j while the mixing strength is defined mainly by 
the overall normalization of distributions to the number 
of excess events observed in the experiments. The anal- 
ysis includes the following constraints: i) the number of 
excess events can vary within the 2a range, ii) the Uh life- 
time has to be less than 2 x 10~ 9 s, iii) the Vh should be 
heavier than 40 MeV to avoid its production in the KAR- 
MEN experiment; and iv) the number of excess events in 
the LSND detector in the energy region 60 < E vis < 150 
MeV with a recoil neutron should be <10 events. The 
latter constrain came from the upper limit on the num- 
ber of events with > recoil neutrons observed by LSND 
in this energy region (38[. 

The summary results are shown in Fig l24l together with 
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FIG. 18: Distributions of the excess events in the MiniBooNE 
detector from the Vh — > jv decay reconstructed as V B CC 
events as a function of E® E for |(7 M h| 2 = 3 x 1CT 3 , m„ h = 40 
MeV, and r Vh — 10~ 9 s, and for different values of the asym- 
metry parameter a. The dots are experimental points for the 
excess events in the MiniBooNE detector. Error bars include 
both statistical and systematic errors 0]. The comparison of 
the distributions with the experimental data yields a \ 2 of 
8.2 (a=-l), 7.1 (a=0), and 6.7 (a=l) for 9 DF. 



the constrains from the it ^2 , K^2 experiments and the 
limit obtained from the recent results of the TWIST ex- 
periment on precision measurements of the Michel spec- 
trum in muon decay, see Sec. VI. The fit results suggest 
that the Vh mass is in the region 40 < m„ h < 80 MeV for 
the Dirac Vh with an asymmetry parameter a = — 1, and 
in the region 40 < m Vh < 70 MeV for the case a = 0. The 
X 2 contribution from the MiniBooNE data is smallest 
for the whole mass range. As expected, for both cases 
the main contributions to x 2 are from the MiniBooNE 
data. For higher Uh masses, preferred by the Mini- 
BooNE data, the region of allowed It^l 2 moves to- 
wards smaller values , while it is cut by the LSND con- 
strain iv). 

The LSND results strongly restrict the allowed Vh mass 
region and exclude solutions with m„ h > 80 MeV, which 
are favored by the MiniBooNE data. The analysis 
gives a 14% \ 2 probability for compatibility between the 
LSND and MiniBooNE data and the Vh — > ~jv interpre- 
tation, demonstrating a reasonable level of agreement. 

As already mentioned in Sec. II, the angular and en- 
ergy distributions of decay photons are sensitive to the 
type of the heavy neutrino. The analysis shows that bet- 
ter fit results can be obtained provided that the zVs pro- 
duced in the LSND and MiniBooNe experiments by the 
muon neutrino decay radiatively as a left-handed Dirac 
neutrino with the asymmetry parameter a = — 1, while 
the i/ft's produced by the muon antineutrinos decay as 
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FIG. 19: Same as FigflHIfor the 70 MeV v h . The comparison 
of the distributions with the experimental data yields a x 2 of 
8.3 (a=-l), 7.5 (a=0), and 6.3 (a=l) for 9 DF. 



a right-handed Dirac neutrino with the asymmetry pa- 
rameter a = +1. The positive sign of the asymmetry 
coefficient is preferred by the analysis of MiniBooNE 77 M 
data, while the negative sign provides a better fit to the 
distributions from LSND and MiniBooNE data. If the 
iv/j's with such exotic properties exist, that would mean 
that the Vh — > 7^ decay is not CP(CPT) conserving 

EH El. 



VI. LIMITS ON |[/ Wl | 2 AND fi tr 

One might reasonably ask if the mixing strength as 
large as the one shown in Fig[M] is consistent with the 
results of previous searches for the MeV Vh 's [45[ . Below 
we discuss the most stringent limits which came from the 
K meson [46| - and muon decays [5(|- [HI], neutrino 
scattering experiments (54[, searches at the CERN LEP 
[55l l56j , and also from cosmology, astrophysics and atmo- 
spheric neutrino experiments [57|-[6l|. Finally, the limits 
on the transition magnetic moment are also discussed. 



A. Limits from K decays 

It is well known that heavy neutrino in the mass range 
< 400 MeV can be effectively probed through the two 
body decays of charged kaons [15| . The K meson, which 
normally decays into a \x and a v^, might instead de- 
cay into a [i and a v^. The experimental signature of 
the presence of the decay K + — s- fxi/^ is a peak in the 
muon energy distribution below the normal one from the 
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FIG. 20: Distribution of the excess events in the MiniBooNE 
from the u h decay reconstructed as V e CC events as a function 
of E vis for E® E > 200 MeV, |f/ Mh | 2 = 3 x 10~ 3 , m„ h = 40 
MeV, and T Vh — 10~ 9 s, and for different values of the asym- 
metry parameter a. The dots are experimental points for the 
excess events in the MiniBooNE detector. Error bars include 
both statistical and systematic errors |8j]. The comparison of 
the distributions with the experimental data yields a \ 2 of 
9.5 (a=-l), 7.5 (a=0), and 6.2 (a=l) for 5 DF. 



ordinary decay at the energy 



E„ 



M 2 



2M 



K 



(21) 



The most stringent current experimental limits on It/^/tl 2 
for the Vh mass region below 400 MeV, are summa- 
rized in Fig. (45(1; for a recent review see [26|, [27j . 
One can see, that the limit for the mass region around 
100 MeV, derived from a search for the Vh at KEK 
49], is |£/ M ft| 2 < 10~ 5 . Surprisingly, the neighboring 
(m Uh ; \Ufj,h\ 2 ) region of parameters (fTij|) . (fTBT) favorable for 
the explanation of the LSND and MiniBooNE results re- 
mains unconstrained. The reason for that is because the 



Vh in the mass range of m v 



40 - 70 MeV is out- 



side of the kinematical limits for tt^ decays and is not 
accessible in decay experiments due to experimen- 
tal resolutions and a high background level. For exam- 
ple, to resolve the muon peak of 234 MeV/c from the 40 
MeV Vh and the main peak of 236 MeV/c, a muon mo- 
mentum resolution better than 1% is required. Another 
reason is related to the K —> /ivy and K + — > /i + tt°v 
decays which produce a continuous background to the 
muon momentum distribution below the main peak and 
essentially constrains the sensitivity of the search for the 
Vh mass range < 100 MeV again due to the require- 
ment of a very high experimental resolution. Let us con- 
sider this in detail. In the most sensitive experiment 
performed at KEK [49j |. degraded K +, s were stopped 



FIG. 21: Same as Fig[20]for the 70 MeV v h . The comparison 
of the distributions with the experimental data yields a \ 2 of 
8.5 (a=-l), 7.2(a=0), and 6.4 (a=l) for 5 DF. 



and decayed in the scintillator target. Charged particles 
from K + decays were momentum analyzed by a magnetic 
spectrometer. To achieve high sensitivity to small sig- 
nals, the main background decay modes, K + — > ^ + tt°v 
and K + /i + ^7, were vetoed by using an almost her- 
metic (92% of 4tt) low-energy threshold (~ 1 MeV) Nal 
calorimeter, surrounding the kaon decay target. The veto 
efficiency for the K + — > {i + ir°v decay mode was quite 
high, better than 99%, thanks to the emission of two 
photons. The decay K + — > ^i + vy was difficult to sup- 
press, and about 30% of photons from this decay mode 
were undetected. The reasons for this are the following: 
i) the low-energy photons are preferably emitted along 
the decay muon momentum direction, so they escape un- 
detected; ii) the photo-nuclear absorption cross section 
is high for photon energies ~ 50 MeV; and iii) there is 
an absorption of decay photons due to the presence of a 
dead material in the vicinity of the target. These effects 
results in a high background level which significantly de- 
creases the sensitivity for the Vh masses below ~ 80 MeV, 
as one can see from Fig. [35] Let us now show that taking 
into account the dominance of the Vh — > "fv decay and 
the relatively short Vh lifetime makes existing experimen- 
tal bounds even weaker. Indeed, in these searches it was 



typically assumed that the Vh is a relatively long-lived 
particle, i.e. — Tzjl. <g; i j where L ~ 50 cm is the typical 

P"h "h 

size of the target region. However, if the decay Vh — > *yv 
is dominant and the heavy neutrino is a short-lived parti- 
cle, the Vh would decay presumably in the vicinity of the 
target. In this case, the decay photon would be vetoed by 
the calorimeter, and the event would be rejected. An es- 
timate shows that for a Vh lifetime, as short as T Uh < 10 -9 
s more than 95% of heavy neutrinos would decay in the 
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FIG. 22: Distribution of the excess events in the MiniBooNE 
detector from the decay reconstructed as V e CC events as 



m„ k = 40 MeV, and r Vh — 10 _a s, and for different values 
of the asymmetry parameter a. The dots are experimental 
points for the excess events in the MiniBooNE detector. Error 
bars include both statistical and systematic errors i- The 
comparison of the distributions with the experimental data 
yields a x 2 of 3.1 (a=-l), 2.7 (a=0), and 3.3 (a=l) for 4 DF. 



targ et region or inside the calorimeter of the experiment 
49], producing a photon with an energy well above the 
veto energy threshold. Because of this self- veto effect, 
the limit \U ^S 1 < (2 — 4) x 10~ 5 for the v h mass around 
~ 80 MeV, could be worsened by more than an order 
of magnitude, and thus, would be in the region ~ 1CP 3 
close to values from (|15[) . Thus, it would be important to 
perform an "open mind" search for heavy neutrino in a 
wider mass range, including the region around 80 MeV. 



B. Muon decay constraints 

If a heavy neutrino with mixing into in the range 
of Eq.tfTof exists, it would notably change the shape of 
the Michel spectrum of the ordinary muon decay, which 
is well predicted in the standard model. This gives the 
possibility of using the results of precision measurements 
of the Michel spectrum in the ordinary muon decay in 
order to probe the possible existence of a heavy neutrino 
fl6| . The relatively free of theoretical uncertainties limit 
on mixing IJT^I 2 for Vh masses in the range from 30 to 
70 MeV was originally set in [13, 49] by using results of 
the p parameter measurement by Derenzo [531 ] . Follow- 
ing Shrock [l6j|, we use the idea that an admixture of 
a heavy neutrino to the Michel spectrum could signifi- 
cantly alter the p parameter, resulting in an effective p 
parameter p e f f that is different from the canonical value 



FIG. 23: The same as in Fig[22] for the 70 MeV v h . The 
comparison of the distributions with the experimental data 
yields a \ 2 of 2.5 (a=-l), 2.7 (a=0), and 2.6 (a=l) for 4 DF. 



p = 0.75. Hence, one can extract limits on the mix- 
ing iC/^ftl 2 from comparison of the measured ( p exp ) and 



< 



effective p e f f values by requiring \p eff - p exp \ ^ a exp , 
where a exp is the error of the measurements. Figure 
1261 shows the dependency of the effective p e f / parame- 
ter as a function of m Uh for the region of interest from 
40 to 80 MeV, for several values of If/^ 2 obtained by 
the fit of the Michel spectrum. The two sigma bands 
from the measurement of the p parameter by Derenzo 
53], p = 0.7518 ± 0.0026, and from the recently re- 
ported pre cision measurements by the TWIST Collab- 
oration [13, p = 0.74977±0.00012(stot)±0.00023(sysi), 
or p — 0.74977 ± 0.00026 with all errors combined in 
quadrature, are also shown. The theoretical expressions 
for the mixing of heavy Dirac neutrinos in muon decay 
can be found in (l6| . and, with radiative corrections in- 
cluded, in 5l| . Our new 2<j limit shown in Fig. [24] is 
derived for the Vh masses in the range from 40 to 80 
MeV by using the results of the TWIST experiment. For 
very large masses the limit is less restrictive because the 
Vh contribution is hi ghly suppressed. Comparison of the 
bounds obtained in |l6l l52j and in the present work by 
using the results of the measurements of the p parameter 
by Derenzo [53j shows good agreement. 



C. Bounds from neutrino scattering experiments 

Next we consider the neutrino experiments that 
searched for Vh decays. The direct searches for the ra- 
diative Vh — > decay were performed for the mass re- 
gion < 1 MeV (45[. Others experiments searched for 
heavy neutrino decaying into charged particles in the fi- 
nal state, such as e.g. v% ~^ eev, jxev, pmv, ... None of 
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FIG. 24: The 2a allowed region (dark areas) in the 
{m Vh ; | C/ M /i | 2 ) parameter space obtained for different values 
of the asymmetry parameter a from the combined analysis of 
LSND and MiniBooNE and 77 M data. The areas excluded 
by the 717,2 and decay experiments [45|. and the exclu- 
sion region obtained in the present work from the results of 
precision measurements of the muon decay parameters by the 
TWIST experiment [Hoi ] are also shown; see Sec. VI. 



these experiments has reported a bound on the mixing 
strength \Ufj,h\ 2 , or on the combination \U^h\ 2 ^tn for the 
radiative Vh — > jv heavy neutrino decay. The experimen- 
tal signature for the Vh decaying into charged particles is 
quite clean. The selection of two tracks originating from a 
common vertex with nonzero invariant mass makes these 
searches almost background-free. In contrast, to search 
for an excess of a single converted photons from the radia- 
tive neutrino decay is more difficult. At high energies the 
background level from the ir° decays and bremsstrahlung 
photons is high. The uniqueness of the LSND and Mini- 
BooNE experiments is that they run at low energies when 
the production of the u^'s is still possible and the back- 
ground level is relatively small due to the high fraction 
of v^NCQE events used for the production of v^s. 

The best limit |[/ m(i | 2 < 1CT 5 - 1(T 6 for the mass re- 
gion m Vh ~ 40 — 80 MeV was derived from a search for 
e + e~ v decays in the PS191 beam dump experiment 



54j . It was assumed that this decay mode is 



at CERN 

dominant and that the Vh is a relatively long-lived parti- 
cle, i.e Lm " h <c 1, where L ~ 1.4 x 10 2 m is the distance 

between the target and the detector. Other decay modes 
with charged particles in the final state, such e.g. as 
Vh — > A* 71 ") MM 17 ) f iei/ are forbidden by the energy conserva- 
tion. The PS191 detector consist of a 12 m long decay 
volume, eight chambers located inside the volume to de- 
tect charged tracks and followed by a calorimeter. The 



FIG. 25: Bounds on the muonic mixing strength \U^h\ 2 of 
the heavy neutrino vs its mass from range measurements 
El, the K + — > fj, + wv decay search experiment [47| . and from 
heavy neutrino searches in 7r — s> \w 48| and K — > [iv [49| de- 
cays. The arrows show the unconstrained LSND-MiniBooNE 
mass window. 



decay volume was essentially an empty region filled with 
helium to reduced the number of ordinary neutrino in- 
teractions down to ~ 100 events, with a total amount of 
dead material around 3.6 g/cm 2 . The events searched for 
in the experiment were requested to consist of two tracks 
originating from a common vertex in the "vacuum" part 
of the Vh decay volume and giving rise to at least one 
shower in the calorimeter. 



Consider now our case, e.g. with \U I 



fih\ 



3 x 10" 



m Vh = 40 MeV and the Vh lifetime T Vh = 10 s. Be- 
cause of the larger mixing the Vh flux from the K decays 
in flight would increase by a factor ~ 10 3 — 10 4 . However, 
several new suppression factors have to be taken into ac- 
count. First, the Vh flux would decrease by a factor ~ 30 
due to the more rapid decay of the v^s. Second, the 
experimental signature for the Vh — > "fv decay would be 
an e + e~ pair from the conversion of the decay photon in 
the decay region. However, the opening angle of the e + e~ 
pair is ~ m e /E e + e - < 10~ 3 rad, which is too small to 
be resolved in the detector, and thus the event would be 
misidentified as a single track event. Such event would be 
rejected. The rejection factor is estimated to be > 10~ 2 . 
Third, the average probability of the photon conversion 
with the vertex located in the low density decay region 
( not in a chamber) is as small as < 10~ 2 . Finally, the 
total number of signal events in PS191 would decrease 
by a factor of more than ~ 10 2 compared to the num- 
ber of events expected for a long-lived Vh 's produced and 
decaying through the mixing IL^/J 2 = 10 -5 . Note, that 
for the above values of \U^h \ 2 , m v h and T Vh , the branch- 
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FIG. 26: Dependence of the parameter p e // as a function of 
the vn,v h mass determined for the decay fi — > ev e fh for several 
values of |t/ M /jj 2 shown in the plot. The 2 a bands around the 
central values of the p parameter measured by Derenzo [H^] 
(dotted lines), and by the TWIST experiment [13] (dashed 
lines) are also shown. 



ing fraction of the direct Vh decay into the e + e~v final 
state is found to be Br(vh —¥ e + e~v) < 10~ 5 26], which 
is also small enough to produce a detectable excess of 
e + e~ events in the PS191 experiment. In this estimate 
the average Uh momentum is < p Vh >~ 4 GeV, the decay 
region length is I — 12 m and the typical photon energy 
is ~ 2 GeV 0. 

D. LEP constraints 

Next we consider bounds from LEP experiments 
pH ]. For the mass region around 50 MeV, the model- 
independent limit from the searches for the Z — > Wh 
decay is \U^h\ 2 J$ 10~ 2 , (see e.g. [55j |) which is compati- 
ble with Eq.(|15p. Direct searches for radiative decays of 
an excited neutrino v* — \ jv produced in Z — > v*v de- 
cays have also been performed [45] |. The best limit from 
ALEPH is H] 

Br(Z -> vv*)Br{v* -> 7 z/) < 2.7 x 10~ 5 . (22) 

As the experimental signatures for the v* — > jv and 
Vh — > jv decays are the same, we will use this bound for 
comparison. The number of expected Vh — ¥ ji> events 
in ALEPH is proportional to Br(Z —> Wh)Br{vh —> 

ju)[l - exp(- m ^ )], with !~lm and p Vh ~ 45 GeV. 



Taking into account ^(z'-Xvt) — \ u ^h\ 2 an d using (J22J, 
we find 

\U^x mVh[M .f ] < 4.8 X 10 9 . (23) 

For the mass range (fT7j)) using Eq. dTBl results in 

Tu h > Kr 11 - 10" 10 s, (24) 
which is consistent with (fl~7)) . 

E. Bounds from cosmology, astrophysics, and the 
Super-K experiment 

Although a detailed analysis of the cosmological and 
astrophysical constrains on the properties of heavy ster- 
ile neutrinos is beyond the scope of this work, let us 
briefly discuss some of them. The most stringent bounds 
|tVl 2 < 10" 10 - 10~ 3 for the v h 's in the MeV mass 
range were obtained from the prim ordial nucleosynthesis 
and SN1987A considerations 57-60]. These limits are 
typically valid under assumption that the Vh is a rela- 
tively long-lived particle with the dominant decay mode 
Vh — > ve + e~ into an active neutrino (y ei v^,v T ) and an 
e + e~ pair. In this case, for the required mass range 40 
- 80 MeV and mixing |C/ M /i| 2 < 10~ 2 , the v^ L lifetime es- 
timated from Eq. pTJl) is r Vh > 10~ 2 s, which is about 7 
orders of magnitude longer compared to the one required 
by (HID. 

Another independent constrain on |E/^fr| can be set 
based on the nonobservation of atmospheric sterile neu- 
trino decays by the Super-K experiment [60(. In this 
work it is assumed that heavy neutrinos could be co- 
piously produced in the Earth's atmosphere and could 
decay inside the Super-K detector, generating an excess 
event rate 60]. The requirement for this rate to not ex- 
ceed the rate of events observed by the experiment results 
in upper limits on mixing strength | C/^ft. | 2 < 10~ 5 — 10~ 4 
for the mass region 40 < m Vh < 80 MeV (see Fig. 6 in 
Ref.[6(|). In these calculations it is assumed that the typ- 
ical Lorentz 7 factor of heavy neutrinos is < 10. Taking 
into account (fT5]) results in a Vh decay length of the or- 
der of / < 10 m. Assuming an average distance between 
the Vh production region and the Super-K detector of 
L ~ 1 km gives a very high i>h flux suppression factor of 
exp(—L/l) -c 10~ 5 . Thus, one can see that the stringent 
bounds from cosmology, astrophysics, and the Super-K 
experiment arc evaded due to the short lifetime of the 
i>h in accordance with (|18p ; for more detailed discussions 
of the bounds on heavy neutrinos from cosmology and 
astrophysics, see e.g. [57l[6l|. 

F. Limits on p tr 

For the light neutrino mass m v << m Vh using Eq.Q 
the Vh lifetime due to a transition moment fj.tr is given 
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by 

l = a ( ^ )9( ^ )a (25) 

The requirement (II 8[) for the Vh 7^ decays to occur 
mostly inside the MiniBooNE fiducial volume results in 

far > 3 x 10~Vb- (26) 

For m Vh ~ 40 - 80 MeV and fi tr > 10~ 8 /is the ra- 
diative decay is dominant, Br{vh — > 71/) > 0.9. Di- 
rect searches for heavy neutrino decays were performed 
by many experiments However, none of these ex- 

periments has reported a bound on the mixing strength 
\Ufih\ 2 or on the product {U^h^fat, for the radiative de- 
cays of heavy neutrino in the mass region 40-80 MeV. 
The mixing |t/ufc| 2 would result in a contribution to the 
effective magnetic moment, [ifff ~ \Ufj,h\ 2 far — (0.4 — 
4.0) x 10 _10 /kb, due to the nonzero Vh magnetic moment. 
This contribution is below the best direct LSND exper- 
imental limit derived from the muon neutrino-electron 
scattering [x e J* < 6.8 x 10~ 10 /is [HI]. However, in this 
particular case the LSND limit is not directly applicable 
to the Vh magnetic moment as the limit was obtained for 
the DAR z7 M , which, as discussed in Sec. II, cannot pro- 
duce Vh in the LSND experiment in v^ scattering, due to 
its heavy mass. 

Consider now again bounds from LEP experiments 
[45| . For the mass region around 50 MeV, the model- 
independent limit from the searches for the Z — > Wh 
decay is |£/ M /i| 2 < 10~ 2 , (see e.g. [55| ) which is com- 
patible with Eq . (fTS")) . Consider the constraint (f2"2"j) from 
direct searches for radiative decays of an excited neu- 
trino v* — > 72/ produced in Z — > v*v decays The 
number of expected v\ x — > 7^ events in ALEPH is propor- 
tional to Br(Z -> vv h )Br(v h -> -yv)[l - exp(- ^ m " h )], 

h Vh 

with |~lm and p Uh ~ 45 GeV. Taking into account 
B Br Z (zZl) ~\U»h\ 2 and using Eq.©, we find 

IM 2 x(^) 2 <1.9xl0- 16 . (27) 

Using Eq. (fl"5"j) results in (j, tr < (2.6 — 1.4) x 10~ 7 [i B , which 
is consistent with Eq. (f26|) . 

The limit on the u t r b etween the Vh and the v^ has 
been obtained in Ref . [631] . based on the idea of the Pri- 
makoff conversion v^Z — > VhZ of the muon neutrino into 
a heavy neutrino in the external Coulomb field of a nu- 
cleus Z, with the subsequent Vh — > jv decay. By using 
the results from the NOMAD experiment [64], a model- 
independent bound [i% r < 4.2 x 10~ 8 /jb was set for the Vh 
masses around 50 MeV (see Table 1 and Fig. 2 in Ref. [63]), 
which is also consistent with Eq. (j26[) . Values of for larger 
than 10~ 8 /i;3 for the m Vh > 40 MeV could be obtained 
e.g. in the framework of the Zee model [22| . 



VII. PROPOSED SEARCHES FOR HEAVY 
NEUTRINO 

In this section we propose experimental searches for 
heavy neutrinos in the decay and muon neutrino 
interactions. The sensitivity of the proposed experiments 
is expected to cover the region of the LSND-MiniBooNE 
parameter space shown in Fig |24l A discussion of the 
possible search for heavy neutrino in muon decays will 
be reported elsewhere. 

A. Search for the Vh in K decays 

As discussed in Sec. VI, the existence of heavy neutri- 
nos with masses < 400 MeV can be effectively probed by 
searching for a peak from the v% in two body K + — > 
decays of charged kaons [l5| . Depending on the experi- 
mental method, one could also search for a peak in the 
missing mass distribution corresponding to the mass of 
the heavy neutrino. The number of K + — > /ivh events 
in the peak is defined by the mixing lE/^l 2 and by the 
phase space and helicity factors which depend on the Vh 
mass [15| . For the mass interval m Vh ~ 40 — 80 MeV the 
chirality-flip is mostly due to the sterile neutrino mass 
which results in 

T(K -+ fivh) « T(K -> ^)\U^\ 2 (^) 2 . (28) 

Using Eq. ([15]). we find that the branching fraction of 
K fiVh is in the experimentally accessible range: 

Br(K ->■ fiv h ) « 10~ 4 - 10~ 3 (29) 

for heavy neutrino masses in the range 40 - 80 MeV. 
There are two advantages to searching for the Vh in the K 
decay peak experiments. First, an observation of a peak 
in the muon and/or the missing mass spectra in K decays 
gives unambiguous evidence for the existence of heavy 
neutrinos. Second, the expected number of signal events 
occurring in a detector, and hence the sensitivity of the 
search, is ex IC/^/J 2 , as it follows from (|2"8")l . In neutrino 
scattering experiments the Vh decay signal rate is either 
proportional to IC^I 2 • \U^h\ 2 , or, if the dominant decay 
is Vh — > jv, proportional to 1/7^ | 2 ■ H 2 r and thus, is more 
suppressed. Here, the first term |/7 M ?i| 2 appears through 
the heavy neutrino production in a target, and the second 
term, (CT^hl or /i 2 r , through the heavy neutrino decay 
in the detector. Note that in our particular case the 
sensitivity of the LSND and MiniBooNE experiments is 
°c I 1 2 because the Vh is a short-lived particle due to 
the large value of fo r , which decays in the detector target 
volume with the probability ~ 1. 

As discussed in Sec. VI, the major physics background 
to the experiments searching for the Vh peak in de- 
cays at rest is the radiative kaon decay K + — > fiv^y which 
has a branching fraction of about 1.5 % for photon en- 
ergy above 10 MeV 45]. This background results in an 



21 



admixture of a continuous spectrum to the muon mo- 
mentum distribution below the main peak and essentially 
constrains the sensitivity of the search for the mass 
range < 100 MeV. To improve the sensitivity, this back- 
ground decay mode has to be suppressed by increasing 
the detection efficiency of the decay photons. Experi- 
mentally, improvement of the photon efficiency for the 
searches with the K decays at rest is difficult due to the 
limitation factors discussed in Sec. VI. 

Here, we propose the use of K decays in flight to im- 
prove the sensitivity against this background source. A 
substantial increase in the detection efficiency of radia- 
tive photons could be obtained by using the K decays 
in flight at high energies. In this case, the vast major- 
ity of decay photons would be within the geometrical 
acceptance of the detector because they are distributed 
within a narrow cone with the maximal photon angle of 
the order 7 ~ nix / Ek — 7 mrad for a kaon energy of 
Ek = 70 GeV, as schematically illustrated in Fig. [27lA 
Thus, the detection efficiency of high energy decay pho- 
tons from the decay K + — > fi + v^ in an electromagnetic 
calorimeter is expected to be almost 100%. If the is 
a relatively long-lived particle, with a lifetime > 10~ 10 s, 
then it would rarely decay in the experiment, as the v>h 
average decay length of > 300 m is much bigger than the 
typical length of a decay volume ~ 100 m. The detection 
of a muon and a photon in the final state would unam- 
biguously signal the detection of the radiative K decay, as 
shown in Fig. I27A . or another background decay mode. 
This would reduce the background significantly and allow 
the measurement of the muon energy distribution with 
a higher sensitivity. If the i>h is a short-lived particle, 
with a lifetime ~ 10~ n [see Eq.(|2"4"|] and a correspond- 
ing decay length of about 30 m, then the detection of a 
muon and a photon in the final state would mean either 
the detection of a background process or the detection 
of the signal from the decay —¥ jf, as shown in Fig. 
[27B. In this case, one could still suppress the background 
by rejection of the fv) observed events at the cost of the 
signal efficiency loss. To avoid this reduction, one could 
try to identify signal events by using the fact that the 
observed photon is originated from a secondary vertex, 
which is displaced from the primary one by a large dis- 
tance, provided that precise measurements of the photon 
directionality can be done. 

A good example of an experiment where the proposed 
search could be performed is the NA-62 at CERN [fjjj. 
The experiment is running at a kaon energy of 74 GeV. 
The detector is well equipped to identify and measure the 
momenta/energy and directions of the charged particles. 
The photons are precisely measured with a LXe electro- 
magnetic calorimeter. To evade the KEK limit for the 
region m Vh < 80 MeV the Vh lifetime is required to be in 
a slightly more restricted range T Vh < 5 x 10~ n s. For 
example, for r Vh ~ 3 x 10 s the most of the i/^'s would 
decay in the K decay volume [65[ , thus producing a veto 
signal in the LXe calorimeter. Other experiments capa- 
ble of searching for the Vh — > 7^ decay with their existing 
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FIG. 27: Schematic illustration of an experiment to search for 
radiative neutrino decay in decays in flight at high ener- 
gies: A) the main background decay K — > fj,wy is suppressed 
because of the high detection efficiency of decay photons in 
the electromagnetic calorimeter (ECAL) due to the Lorentz 
bust; B) if the Vh is a short-lived particle, a part of the pho- 
tons from the Vh — > 7^ decay is also detected. See text. 



data are the E787 and its upgrade, the E949, at BNL [Hj], 
or ISTAR+ at IHEP \§7§. The former is equipped with 
an almost 4-7T veto electromagnetic calorimeter, allowing 
good rejection of photons from background K + — > [i + ir b> 
and K + — > [i + h>j decays. 



B. Search for the decay v% — \ 'yv in NC neutrino 
interactions 

As discussed above, in order to search for an excess 
of single converted photons from the radiative neutrino 
decay Vh — > 7^ in high energy NC neutrino interac- 
tions the background, mainly from the decays of 7r°'s and 
bremsstrahlung photons produced either in the primary 
vertex or in the secondary particles interactions has to be 
eliminated. To suppress the background and to detect a 
clean and convincing sample of converted decay photons 
one can perform a neutrino "beam dump" experiment the 
main idea of which is illustrated in Fig. 1281 A neutrino 
detector is subdivided into two parts. The first part is 
an active absorber part, and the second one is the decay 
region for the detection of converted photons from the 
decay Vh — > iv. The secondary particles from z/^NC in- 
teractions in the detector are absorbed in its first part. 
Heavy neutrinos produced through the muonic mixing 
penetrate the dump and decay into a photon and a light 
neutrino in the second downstream part of the detector, 
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with a subsequent photon conversion. The experimen- 
tal signature of the decay — > 7^ is the appearance of 
a single e + e~ pair originating from a secondary vertex 
displaced from the primary one at a distance L signifi- 
cantly larger than the detector nuclear interaction length, 
L 3> Ao. The main background sources for this setup are 
expected from the secondary neutrons and/or K^s pen- 
etrating the dump and producing 7r°'s either in hadronic 
secondary interactions or in decays in flight in the target. 
The decays of these 7r°'s could be misidentified as a sin- 
gle decay photon event. The suppression of these back- 
grounds can be achived by increasing the n, absorp- 
tion in the first part of the detector simply by increasing 
its length or by selecting events with larger L. Obviously, 
the precise identification of the electromagnetic nature of 
the signal event is of great importance for this search. In- 
terestingly, if the event excess is indeed originated from 
the converted photons, the proposed search could also 
distinguish whether the excess events are produced by 
photons from the i>h — » decays or by ones emit- 
ted from the primary vertex due to anomaly- mediated 
neutrino-photon coupling, as discussed in [68[. The (al- 
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FIG. 28: Schematic illustration of the proposed neutrino 
experiment to search for radiative neutrino decay in u^NC 
interactions. The electromagnetic and hadronic secondaries 
from the v^NC event are absorbed in the initial part of a neu- 
trino detector serving as a dump. Heavy neutrino produced 
through muonic mixing (see Fig[T} penetrates the dump and 
decay into a photon and a light neutrino in the downstream 
decay region of the detector. The experimental signature of 
the Vh — > decay is the appearance of a single high energy 
e + e~ pair from the conversion of the decay photon at a dis- 
tance L from the primary vertex significantly larger than the 
detector nuclear interaction length L 2> Ao. Precise identi- 
fication of the electromagnetic nature of the excess events is 
crucial for this experiment. 

most) ideal detector to search for the decay Vh — > 7^ is 
a detector similar to the NOMAD one The NO- 

MAD is equipped with a forward calorimeter (FCAL), 
where the secondaries from high energy neutrino inter- 
actions in FCAL could be dumped. In addition, it has 
the excellent capability of identifying and reconstructing 
converted photons due to the low mass target located in 
a magnetic field. An example of the reconstruction of 
two conversion e + e _ pairs from the decay tt° — > 2j can 
be found in Ref. [70] . One disadvantage of the detector is 
the short length of the tracking part. The overall detec- 



tion efficiency of the production, decay in flight with 
the subsequent conversion of the decay photons into an 
e + e _ pair and the reconstruction of the conversion pairs 
is expected to be low. The advantage of NOMAD is its 
great capability to measure the e + e~ pair directionality 
with a precision of (I — cosQ e + e -) < 10~ 5 71]. This 
will allow an effective suppression of converted photons 
originating from the primary vertex; see the discussion 
in Ref. [7l| on tt° reconstruction. 

Another experiment capable of searching for the de- 
cay Vh — > 7^ is the ICARUS T600, which is currently 
taking data at the CERN-Gran Sasso neutrino beam 
[72l ]. The detector is composed of two identical adjacent 
T300 half- modules filled with liquid argon (LAr). A de- 
tailed description of the apparatus can be found in [72j . 
Each T300 half-module has the following internal dimen- 
sions: 3.6 x 3.9 x 19.9(length) m 3 . LAr has a radiation 
length of Xq = 14 cm and a nuclear interaction length 
of A = 83.6 cm, and therefore provide good electromag- 
netic and hadronic secondary absorption and detection 
capabilities for the proposed search, assuming that the 
length of the decay region is L > 10 m^ Xq, Xq . 

The number of Vh — > events in ICARUS can be 
estimated as follows 

AN Vh ^ v ~ N NC \U^ h \ 2 P dec P abs e (30) 

where Nnc — 10 3 is the number of the detected neutral- 
current events, and Pd ec {— 0.4),P a f, s (~ l),ande(~ 0.7) 
are the probabilities for the Vh decay in the detector fidu- 
cial volume and decay photon conversion, and the overall 
detection efficiency of the e + e~ pair, respectively. In this 
estimate the average i>h momentum is < p Vh >~ 10 GeV, 
T„ h < 10 ~ 9 , and the length of the decay region is L = 12 
m. Finally, we find that the number of expected Vh — > 7^ 
signal events in ICARUS is 

AN„ h ^ JV ~ 6 x 10 2 x \U, lh \ 2 (31) 

For the allowed mixing (see Fig. [24]) , this results in 
AN^-^^v ~ 1 — 3 events. If no candidates are seen above 
the expected background level, ICARUS could set a limit 
on the mixing strength of the order < 10~ 3 , which 

is competitive for the mass range 40-80 MeV with the 
bounds obtained from the TWIST experiment, possibly 
allowing us to rule out the LSND-MiniBooNE parame- 
ter region. Note that the search for an excess of the Vh 
decay events can also be performed in the recently pro- 
posed ICARUS-like experiment at CERN PS [zl, or at 
FNAL with neutrino detectors such as MicroBooNE [74j , 
HiResM;/ [75j], and B00NE (a MiniBooNE near detector) 

VIII. SUMMARY 

In summary, we reexamined neutrino oscillation re- 
sults from the accelerator experiments LSND, KARMEN 
and MiniBooNE. We showed that the LSND evidence for 
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— > 77 e oscillations, its long-standing disagreement with 
the results from KARMEN, and the anomalous event ex- 
cess observed by MiniBooNE in and data, can all be 
explained by the production and decay of a heavy neutral 
lepton. The shape of the excess events in several kine- 
matic variables in the LSND and MiniBooNE and 77^ 
data is found to be consistent with the distributions ob- 
tained within this interpretation, assuming that the v^s 
are created by mixing in neutral-current interactions 
and decay radiatively into 72/. Therefore, our main pre- 
diction is that the excess of events observed in the LSND 
and MiniBooNE experiments originates from the Comp- 
ton scattering or e + e~ conversion of the decay photons 
in these detectors. In this context, the confirmation of 
the photon origin of the excess events by measurements 
with a detector able to distinguish electrons and photons 
becomes a crucial test for this scenario. 

A combined analysis of the energy and angular distri- 
butions of the excess events observed in the LSND and 
MiniBooNE experiments suggests that the Vh mass is in 
the range from 40 to 80 MeV, the mixing strength is 
I £7^ 1 2 ~ 10~ 3 - 10~ 2 and the lifetime r Vh < 10~ 9 s. Sur- 
prisingly, this LSND-MiniBooNE favorable parameters 
window is found to be unconstrained by the results from 
the most sensitive K^, neutrino scattering, and LEP 
experiments. Because of the short Vh lifetime, the con- 
straints coming from cosmological and astrophysical con- 
siderations, as well as the bounds from the atmospheric 
neutrino measurements, are also evaded. We set new 
limits on the mixing |C/^ft| for heavy neutrino masses 
in the range 40 to 80 MeV by using results on precision 
measurements of the Michel spectrum by the TWIST ex- 
periment. We also discuss the most natural model for the 
Vh — > decay through the transition magnetic moment 



between the Vh and the light neutrino and show that the 
obtained values |f/ M ^| 2 ~ 10~ 3 - 10~ 2 and pL tr > 10~Vb 
do not violate bounds from previous experiments. 

The results obtained provide a strong motivation for 
a sensitive search for the Vh in a near future K decay 
or neutrino experiments. We propose such experiments 
with the expected sensitivity to cover the region of the 
LSND-MiniBooNE parameter space and notice they fit 
well in the existing and planned experimental programs 
at CERN or FNAL. The radiative heavy neutrino decay 
could be present in various extensions of the standard 
model and, thus, could enhance the reported motivations 
to search for this process. We note that an analysis of 
the excess of events due to the v% — > jv decay may also 
be possible with existing neutrino data; e.g. new results 
could be obtained from NOMAD [II]. 

The reported analysis gives the estimated values of the 
parameters m Vh , IL^/J 2 , and r^h and may be improved by 
more accurate and detailed simulations of the LSND and 
MiniBooNE detectors. It would also be interesting and 
important to have general analysis of the production of 
heavy neutrinos of Dirac or Majorana type, e.g. in v^NC 
interactions, for arbitrary weak couplings, including the 
leptonic mixing and helicity effects. 
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